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ABSTRACT 

Aims. We study the emergence of a non-twisted flux tube from the solar interior into the solar atmosphere. We investigate whether the 
length of the buoyant part of the flux tube (i.e. A) affects the emergence of the field and the dynamics of the evolving magnetic flux 
system. 

Methods. We perform three-dimensional (3D), time-dependent, resistive, compressible MHD simulations using the Lare3D code. 
Results. We find that there are considerable differences in the dynamics of the emergence of a magnetic flux tube when A is varied. 
In the solar interior, for larger values of A, the rising magnetic field emerges faster and expands more due to its lower magnetic 
tension. As a result, its field strength decreases and its emergence above the photosphere occurs later than in the smaller A case. 
However, in both cases, the emerging field at the photosphere becomes unstable in two places, forming two magnetic bipoles that 
interact dynamically during the evolution of the system. Most of the dynamic phenomena occur at the current layer, which is formed 
at the interface between the interacting bipoles. We find the formation and ejection of plasmoids, the onset of successive jets from the 
interface, and the impulsive heating of the plasma in the solar atmosphere. We discuss the triggering mechanism of the jets and the 
atmospheric response to the emergence of magnetic flux in the two cases. 
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1. Introduction 


Active region formation and dynamic phenomena on the Sun 
have been associated with magnetic flux emergence from both 
a theoretical and observational perspective (e.g. |Parker| 1955t 
Tsinganos |1980[ [Riest|| 19821 [Zwaan 1985[ Shibata 1999[ Ar- 


chontis et al.|20(j4[|Schrijver|2()09| l. In numerical flux emergence 
models, a typical initial condition for the subphotospheric mag¬ 
netic held is a twisted magnetic flux tube, which undergoes a 
buoyant rise towards the solar surface. 


The buoyant rise is initiated by removing plasma den¬ 
sity along the flux tube axis in a Gaussian manner (Ap cx 
sxp(-y^/A^)). The length of the buoyant part of the flux tube 
(commonly denoted by d) is an important parameter that affects 
the emergence of the flux tube and the following dynamics. For 
instance, |Fan| ( |200r| ) performed an experiment of a flux tube with 
d = 20 (i.e. ss 3 Mm), which emerged into the solar atmosphere. 
[Manchester et al.| ( |2004| l, used the same initial conditions as [Fan j 
( 2001| l, but with half the length of the buoyant section of the flux 
tube (d = 10). They found that the decrease in d results in more 
efficient plasma draining. In addition, they found the formation 
of a new flux rope at the low atmosphere due to the shearing and 
reconnection of the fleldlines across the polarity inversion fine 
(PIF). [Archontis & Hood| ( |2010| ) studied the effect of different 
field strengths, twist, and d to the formation of magnetic tails. 
For a highly twisted flux tube with high held strength, located at 
the upper convection zone, they changed the buoyant part of the 
flux tube from d=10 to d=20. This resulted in less magnetic ten¬ 


sion on held lines of the higher d during the initial emergence, 
leading to relatively faster emergence. 

Another parameter that plays a major role in the emergence 
and the dynamics of flux tubes is the twist of the magnetic 
held. In many twisted flux tube models (e.g. Fan 2001 1 Mag- 


jara & Fongcope||2003[ [Archontis et al.j 20051, the held fines 
wind about the axis of the tube in a uniform manner. A com¬ 
mon setup for such axisymmetric flux tubes is to assume that the 
flux tube has a component of its magnetic held along its axis 
(Ba.xis) and a poloidal component that is proportional to Baxis 
(i.e. oc aBaxLs). This parameter a measures the twist of the 
magnetic held per unit of length. It also defines the pitch an¬ 
gle between the magnetic held vector and the axis of the flux 
tube. Highly twisted flux tubes with pitch angle greater than 45° 
would become kink unstable ( [Finton et al.[1996|[Fan et al.|1999| l. 

Concerning the ability of low-twist magnetic fields to emerge 
above the photosphere, simulations have indicated that the flux 
tubes can become fragmented into two counter-rotating vortices 
during their emergence ( [Schuessf^[1979 1 . This fragmentation 
can be suppressed by the addition of a sufficient amount of twist 
(e.g. Moreno-Insertis & Emonet 1996| l. On tfie otfier hand, by 
performing simulations with different perturbations of specific 
entropy in tfie flux tube, [Abbett et al.[ ( [2000| l found that even 
when the twist of the initial tube is not sufficient, the tube can 
remain coherent provided it has sufficiently large curvature at its 
rising apex. Abbett et'ar] ( 200l| found that the Coriolis force also 
suppresses the degree of fragmentation of the apex of the emerg¬ 
ing flux tube. Also, if the magnetic held is much stronger than 
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(a) Initial Stratification (b) EFR center 




(c) Bipole center 



Fig. 1: (a): Initial stratification of the atmosphere (temperature (T), density (p), magnetic pressure (Pm), and gas pressure (Pg)). (b): 
Instability term at the center of the EFR and (c): at the center of the bipole at ^ = 4.8 Mm. 


the equipartition value of the magnetic field in the convection 
zone, the buoyancy force dominates over the convective down- 
drafts, and the flux tube can emerge unimpeded ( |Fan et al.|2003| l. 
Failed emergence of flux tubes with low twist have been reported 
in several studies. Murray et al. ( 2006| l performed experiments of 
flux tubes with A — 20, for different values of twist (a = 0.1 -0.4) 
and found that the highly twisted flux tube (a - 0.4) emerged 
into the atmosphere forming a bipolar region on the photosphere. 
Flux tubes with intermediate twist (a = 0.2) emerged form¬ 
ing two bipolar regions (like the “sea serpent” configuration) on 
the photosphere. Weakly twisted flux tubes (a - 0.1) failed to 
emerge above the photosphere. Instead, they continued to ex¬ 
pand horizontally inside the solar interior. |Toriumi & Yokoyama| 
( 201 l| l also studied the emergence of flux tubes with A - 20 
(and similar initial conditions to [Murray et al.|2006|l for a wider 


range of twists (a — 0.05 - 0.5) and found similar results to Mur- 
ray et al. ( 2006| l. Those results indicate that weakly twisted flux 
tubes struggle to emerge above the photosphere. 

The effect of A (which affects the curvature of the emerging 
part of flux tube) on the emergence of weakly twisted flux tubes 
was studied by [Archontis et al. (|2013| ). They used similar ini¬ 
tial conditions to Murray et al. ( 2006|) and Toriumi & Yokoyama 


(201 l|l for a flux tube with a - 0 . 1 , but assumed half the length 


of the buoyant part of the flux tube (d = 10). In this case, the 
weakly twisted flux tube did not fail to emerge above the photo- 
phere. Instead, the flux tube formed two bipolar regions on the 
photosphere (similar to the a = 0.2 case of [Murray et al.|2006| l. 
Above the two bipoles, the expanding magnetic held formed two 
magnetic lobes. The interaction of the lobes triggered several re¬ 
connection jets and created an envelope held that confined them. 
The shearing motion and reconnection of held lines along the 
PIL led to the formation of a highly twisted flux rope, which 
remained at the low atmosphere during the running time of the 
simulation. 


In this paper, we advance the work by [Archontis et al.| ( |2013| l 
on the emergence of weakly twisted fields by adding two key el¬ 
ements: Firstly, we use a non-twisted flux tube and secondly, we 
study the effect of A on the overall dynamics of the emerging flux 
system. The use of a non-twisted flux tube (a - 0) is chosen be¬ 
cause non-twisted fields are more susceptible to distortion and/or 
fragmentation and, consequently, to an overall failed emergence. 
Furthermore, we want to examine whether the emergence of a 
non-twisted magnetic held can create twisted magnetic struc¬ 
tures in the solar atmosphere, similar to the post-emergence flux 
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ropes formed at the weakly twisted case by [Archontis et al.[ 
( |2UT3] |. 

This paper is organized in the following manner. In Section]^ 
we describe the initial conditions of the experiments. In Section 
1 ^ we study the case of T = 10 in terms of its initial emergence, 
the dynamics of the interaction between the formed magnetic 
bipoles, the coronal response, and the energy and flux transfer 
into the atmosphere. In Sectionwe study the d = 5 case with 
respect to the same aspects and compare the results with the A = 
10 case. In Section]^ we summarize the results. 


2. Initial conditions and numerical code 


In order to study the emergence and evolution of the non-twisted 
flux tube, we solve the 3D time-dependent, resistive, compress¬ 
ible MHD equations in Cartesian geometry. We use the Lare3D 
code of |Arber et al.[ ( [200T| l to solve numerically the equations. 
The equations in dimensionless form are 


— H- V ■ (pv) = 0, 
of 

= -V ■ (pvv) -H (V X B) X B - VP -H pg -H V ■ S, 
ot 

= -V ■ (pev) - PV ■ V H- Qjoule + Gvisc, 

— = V X (v X B) -I- nV^B, 

dt V 7 , 

P 

(r-i)p’ 


( 1 ) 

( 2 ) 

(3) 

(4) 

(5) 


where p, v, B, and P are density, velocity vector, magnetic held 
vector, and gas pressure. Gravity is included. We use constant 
explicit resistivity ( 77 ) such that the magnetic Reynolds number 
is Rm-IQO. We assume perfect gas specific heat of y = 5/3. Vis¬ 
cous heating Gvisc and Joule dissipation Gjouie are also included. 

The dimension constants we use are, density pc = 1.67 x 
10'^ g cm“^, temperature Pc = 5100 K, pressure Pc = 7.16 x 
10 ^ erg cm“^, length = 180 km and magnetic held strength 
Pc = 300 G. From the above units we get velocity vq = 
2.1 km s“' and time to = 85.7 s 

The numerical box consists of an adiabatically stratified sub- 
photospheric layer at -3.6 Mm < z < 0 Mm, an isothermal 
photospheric-chromospheric layer at 0 Mm < z < 1.9 Mm, a 
transition region at 1.9 Mm < z < 3.2 Mm and an isothermal 
coronal at 3.2 Mm < z < 25.2 Mm. 



























































P. Syntelis et al.: Emergence of non-twisted magnetic fields in the Sun 


(a) 



(b) 



Fig. 2: (a): Magnetogram at the mid-photosphere (z — 0.8 Mm) at f = 169 min. For simplicity, we call A,D the outer polarities and 
B,C the inner polarities, (b): Field lines traced from the bipoles at f = 169 min. 


We assume to have a field-free atmosphere in hydrostatic 
equilibrium. The initial distribution of temperature (T), density 
(p), gas (Pg) and magnetic (Pm) pressure is shown in Fig. 

For the initial magnetic field, we assume a horizontal cylin¬ 
drical non-twisted magnetic flux tube located at 2.1 Mm below 
the photosphere. The axis of the flux tube is oriented along the 
^-direction, so the transverse direction is along x and height is in 
the z-direction. 

The axial field of the tube is 


By = Boexp(-r^/P^), 


( 6 ) 


where R - 450 km is a measure of the tube’s radius and r the 
radial distance from the tube axis. In order to make the flux tube 
to emerge, we assume a density distribution similar to the work 
by Archontis et al. (20041, which makes the middle part of the 
flux tube less dense and thus buoyant. The density deficit is 


Ap = ^^p(z)exp(-i/2/T2), 

P(z) 


(7) 


where p is the external pressure and pt is the pressure within 
the flux tube (magnetic plus gas). The parameter A is the length 
scale of the buoyant part of the flux tube. In our numerical ex¬ 
periments, we use /I = 5 (0.9 Mm) and 4 = 10 (1.8 Mm). The 
field strength of the flux tube is Bq - 2780 G. 

The computational box has a size of 28.8 x 28.8 x 28.8 Mm 
in the x,y,z directions, in a 457 x457 x457 grid. The resolution is 
such that, in the sites with steep field gradients (e.g. current lay¬ 
ers), the magnetic Reynolds number is (on average) larger than 
the effective Reynolds number. This indicates that the explicit re¬ 
sistivity dominates numerical effects. We assume open boundary 
conditions in the x direction, periodic boundary conditions y di¬ 
rection, and a wave damping zone for z > 26 Mm, which decays 
exponentially the velocity, gas pressure and density fluctuations. 


3. Case of T = 10 

In this section we describe the 4=10 case. We study the emer¬ 
gence of the flux tube through the highly stratified atmosphere 
and the onset of the associated dynamic phenomena in the EFR. 



Fig. 3: Grayscale: B^ component of the magnetic field at the mid¬ 
photosphere (z = 0.9 Mm), ranging from -40 G to 40 G. Colored 
scale: density at i^z-midplane, ranging from 1.7x10“^"^ g cm~^ 
(red) to 1.7x10 g cm“^ (blue). The lines correspond to field¬ 
lines connecting the bipoles (blue), fieldlines of the arcade (red) 
and fieldlines of the envelope field (orange). Time is f = 250 min. 


3.1. Initial emergence of the flux tube 

As the magnetic field emerges from the solar interior, the flux 
tube reaches the photosphere where it becomes compressed. In 
order for the flux tube to emerge above the 
to satisfy the buoyancy instability criterion 
[chontis et al.|2004| ) 

-//p|(logB)>-|^d + k2(l + |j, (8) 

where Hp is the pressure scale-height, B is the magnitude of the 
magnetic field, y - 5/3 is the specific heats ratio,/! is the plasma 
parameter, 6 - -0.4 is the superadiabatic excess and k are the 
perturbation wavenumbers, where k\\ and are the horizontal 
components parallel and vertical to the magnetic field, and 
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photosphere, it has 
(Acheson|1979|[S7 
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Fig. 4: (a): Temporal evolution of maximum (black) and minimum (red) above the photosphere. Vertical dashed lines correspond 
to the emission of the bi-directional jets at f = 242,248,257 min. (b): Maximum temperature (black) and density (blue) around the 
locations of maximum V^. 


is the vertical component. To calculate these wavenumbers, we 
measure the size of the flux tube sections below the photosphere. 
We hnd « 0.5 Mm and /y ^ 1.4 Mm and approximately the 
radius of the tube as the vertical width 4 ~ 0.45 Mm. These es¬ 
timates are taken at f = 154 min. The corresponding wavevector 
components are k\\ ^ 4.5 x 10“* cm“', 1.2 x 10“^ cm“' and 

kj 1.4 X 10"^ cm“4 Since the wavevector term in the insta¬ 
bility criterion is small (4.6 x 10“^^ cm“^) in comparison to the 
jS term (which has minimum values of 10“'^), we neglect it, and 
after dividing the inequality with its right part, we get approxi¬ 
mately: -//p^(logZ?)/(-|j8d) > 1. We refer to the left side of 
this equation as the “instability term”. 


At the center (around .r = 0, = 0, z = 0) of the EFR, 

the field magnitude decreases as the flux tube expands. At the 
areas sideways of the center, along the j/-direction, the held is 
compressed, leading to a smaller gradient of the magnetic held 
and a lower jS. This makes the instability term greater than 1 in 
these areas, while it remains below unity around the center of 
the EFR. Thus, the flux tube does not emerge at the center of 
the EFR, but it forms two side bipoles, as shown in Fig. |^. In 
Fig#>’[5 we plot the temporal evolution of the instability term 
at two positions in the numerical box. In Fig. [^, we calculate 
the instability term at the center of the EFR. We confirm that 
the held does not emerge at this region (Fig. [^). The flux tube 
becomes unstable at two areas along the flux tube axis, at (x, y) - 
(0, +4.8 Mm). There, the instability term becomes greater than 1 
(Fig- calculated aty- 4.8 Mm). This leads to the emergence 
of two bipolar regions. 


At f = 169 min (Fig.|^), the two bipoles have emerged into 
the photosphere. Hereafter, we refer to polarities A,D as outer 
polarities, and C,B as inner polarities. During the evolution, the 
polarities at each bipoles expand along the y-direction. At the 
photosphere, the fleldlines that join the two polarities of each 
bipole are almost horizontal due to the absence of twist. If the 
flux tube was twisted, we would have detected magnetic “tails” 
on the two sides of the PIE in each bipole (e.g. Lopez Fuentes 
|et al.||2000| [Chandra et al.|[2009| ). As the magnetic held rises 
above the photosphere, the magnetized volume of each bipole 
expands and the attached magnetic held lines form two magnetic 
bipoles (Fig.|^). 
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3.2. Interaction of the magnetic bipoles 

In the following, we study the effect of the interaction of the 
magnetic bipoles on the solar atmosphere. Firstly, we focus on 
the formation of plasmoids and jets, and secondly on the transfer 
of plasma and flux into the solar atmosphere. 


3.2.1. Jets and plasma transfer 


After f = 171 min, the two magnetic bipoles come into con¬ 
tact and a current layer is formed at their interface. Eventually, 
the long and thin current layer becomes subject to the tearing 
mode instability. Thus, plasmoids are formed within the inter¬ 
face and later on they are ejected out of it. This induces inflows 
towards the interface and initiates effective reconnection of the 
oppositely directed fleldlines of the bipoles into contact. The dy¬ 
namics of plasmoids in a similar configuration (i.e. interaction 
between two individual twisted flux tubes) was studied by |Ar-| 
chontis et al. ( 2007| l. Eig. shows the overall topology of the 
fleldlines at a later stage of the evolution, when two new mag¬ 
netic flux systems (apart from the two magnetic bipoles) have 
been developed. More precisely, reconnected fleldlines join the 
two outer polarities forming the envelope held (Eig. orange 
lines) and reconnected fleldlines join the inner polarities forming 
an arcade-like structure (Eig. red lines). Because of the emer¬ 
gence of the held and the gravitational draining of the plasma, 
low-atmosphere material accumulates at the interface. Ailing the 
area with cool and dense plasma (Eig. green color around in¬ 
terface in the j/z-plane). 

In Eig. we plot the temporal evolution of maximum 
(black) and minimum (red) V4 above the photosphere. Up to f ^ 
240 min we And moderate upflows/downflows (5-10 km s“') 
due to the emergence of the two bipoles and the draining of the 
plasma at the periphery of the magnetic bipoles respectively. The 
first major upflow occurs at f = 242 min (marked with a vertical 
dashed line). This is associated with the ejection of a plasmoid 
from the interface current layer into the corona. The ejection is 
driven by the tension of the reconnected fleldlines at the inter¬ 
face. The ejected plasmoid is cool 10^A") and dense (1-2 or¬ 
ders of magnitude heavier than the background plasma in the 
corona). During the evolution of the system, several plasmoids 
are formed and eventually being ejected from the 3D interface. 
The ejection of the plasmoids is followed by (fast) reconnection 
of fleldlines at the interface and the emission of a series of bi- 
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Fig. 5; The component of the velocity field on the t/z-midplane. Field lines show the projection of the full magnetic field vector 
onto that plane, (a): Prior to plasmoid ejection at f = 242 min. (b): Jet ejection at f = 248 min. 


directional flows (jets). The occurence of the bi-directional jets 
is shown in Fig. for f > 242 min. 

Fig .|^ shows the maximum temperature and density around 
the locations of maximum V^. We find that the temperature evo¬ 
lution of the temperature is also highly impulsive. In fact, there 
is a very good temporal correlation between the bi-directional 
flows and the sudden temperature enhancements, which is an in¬ 
dication of a reconnection-driven plasma acceleration at the in¬ 
terface. The impulsive behavior shows that there are many jets, 
which are emitted from the interface, during the interaction of 
the magnetic bipoles. The first three high-speed jets occur in the 
period t 240 - 260 min. They have temperatures in the range 
1-1.5 MK and densities between 10“'® - 10“*® g cm“® (Fig.[^) 
and, thus, they might account for soft X-ray jets. The following 
jets (for f > 260 min) are (in general) less energetic, with ve¬ 
locities up to 50 km s“* (with an exception, at f ^ 295 min, 
where V(. ^ 100 km s“*) and temperatures lower than 1 MK. 
The general decay of the onset of high-speed jets is due to the 
gradual decrease in the available flux and energy of the system 
that drives the ejections. The available energy comes only from 
the emerging magnetic field. Thus, this energy is limited and it 
unavoidably becomes exhausted. 

Fig .|^ is a close-up of the interface showing the upward mo¬ 
tion of the plasmoid before the onset of the first bi-directional jet. 
The colorscale shows the vertical component of the velocity field 
and the overplotted fieldlines (black lines) show the projected 
magnetic field vector onto the plane. The fieldlines reveal the 
existence of the plasmoid (magnetic “island” in the 2D plane) at 
the interface. The upflow (red color) illustrates the upward mo¬ 
tion (ejection) of the plasmoid at f = 242 min. The topology 


of the fieldlines show that the ejection is indeed driven by the 
tension of the V-shaped reconnected fieldlines underneath the 
plasmoid. The downward product of reconnection is the emis¬ 
sion of a weaker downflow (blue color) and the formation of the 
O-shaped arcade. 


Fig. shows the emission of a strong bi-directional jet af¬ 
ter the expulsion of the plasmoid out of the interface. Now, the 
fieldlines from the interacting bipoles come into direct contact 
(since the plasmoid in between them has been disposed) at the 
interface and reconnect. This is reminiscent of the plasmoid- 


induced reconnection scenario, suggested by Yokoyama & Shi- 
|bata| ( |1994| l. The V-like reconnected fieldlines, which are re¬ 
leased upwards, accelerate the reconnection outflow (jet). The 
downward-released fieldlines are added to the arcade, which 
thus, undergoes an apparent rising motion and it grows in size. 


Fig. l^-c show the plasma dynamics during the emission of 
the reconnection jet at f = 248 min. The bi-directional flow starts 
at a coronal height (z ~ 4.5Mm). The upward reconnection flow 
runs with a velocity of 190 km s“* and it carries relatively dense 
and hot (> 1.0 MK) plasma into the corona (6a, 6b). It reaches a 
height of ss 10 Mm, where it collides with the envelope field and 
it is diverted into two side flows, which run with velocities of 
ss 68 km s“*. Figure.]^ is a close-up of the reconnection region 
showing clearly the inflows towards the reconnection site, which 
are enhanced after the ejection of the plasmoid, and the profound 
hot, bi-directional jet driven by reconnection. 
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Fig. 6: Density (a) and temperature (b, c) distribution at the vertical i/z-midplane. Overplotted is the projection of the velocity field 
vector onto the plane. Bottom: A close-up of the reconnection site. 


3.2.2. Undulation of the PIL, fieldline topology and structure 
of the jets 

One of the interesting features in this experiment is the evolution 
of the arcade structure at the interface. To study this, we firstly 
show the temporal evolution of the distribution at the middle 
of the photosphere (Fig. Initially (Fig.[^, t = 243 min), 

the inner polarities collide and the PIL in between them be¬ 
comes almost perfectly aligned with the jr-axis at the center of 
the emerging flux region. However, later we find that the PIL 
adopts an undulating shape (Fig.jTjr, f = 286 min). The overplot¬ 
ted velocity field on this plane shows that the deformation of the 
PIL is associated with plasma motions along the y-axis, towards 
the PIL. A similar deformation was reported in the work by |Ar-| 
jchontis & Hood| ( |2010| l: the tails of the magnetic field on the two 
sides of the PIL suffered by converging flows (inflows), which 
were produced due to the rapid expansion of the field above the 
photosphere. Here, we find that a similar process is at work. The 
finger-like segments, which are developed across the PIL, are 
sites with low density and pressure. At this height (z - 0.9 Mm), 
the magnetic field is also weak and the kinetic energy density is 
such that the magnetic field is advected towards the PIL. The low 
pressure at various points is produced due to the rapid ejection 
of several dense plasmoids from along the interface. The plasma 
then moves towards the low-pressure areas (inflows), deforming 
the PIL. 

Figure]^ shows the 3D topology of the fieldlines around the 
arcade before the deformation (f = 243 min). Overplotted is the 
full magnetic field vector at the photosphere, colored by the ver¬ 
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tical (Bj) component (red is positive, blue is negative). The ar¬ 
cade consists of fieldlines that join the inner polarities all along 
the PIL. The fieldines have been traced from random positions 
within the two inner polarities. We note that the fieldlines of the 
arcade towards its two ends (along the jr-direction) are more in¬ 
clined than the fieldlines at the middle of the arcade. This is due 
to the fan-like shape of the interacting fieldlines of the two mag¬ 
netic bipoles. We also note the cusp-like shape at the top of the 
arcade’s fieldlines. This indicates the occurrence of reconnection 
above them. Left and right from the arcade, we find expanding 
fieldlines that belong to the interacting bipoles. Envelope field¬ 
lines are not shown in this figure. 

Figure]^ shows the topology of the fieldlines around the de¬ 
formed PIL at f = 286 min. We note that after the deformation, 
the magnetic fieldlines becomes highly twisted all along the PIL. 
Some segments of the PIL are now oriented along the y-axis. For 
instance, the segment which is outlined by the red rectangular in¬ 
set in Fig.|^, is the same with the segment around x = -0.5 Mm 
and between y = -1 Mm and y - 1 Mm in Fig. |^. Now, the 
relative motion of the inflows on the two sides of this segment 
induce a shearing (Fig. [Tjr, motions towards negative t/-axis and 
positive i/-axis around x = -0.3 Mm, as the positive (negative) 
inner polarity moves towards the negative (positive) polarity). 
The sheared and oppositely directed fieldlines across the PIL’s 
segment reconnect, producing a small twisted flux rope and a re¬ 
connection jet. This process occurs at several (similar) segments 
along the PIL, providing an additional source of twist for the 
fieldlines around the PIL and triggering the onset of more jets 
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Fig. 7: (a): component of the magnetic field at the mid-photosphere {z - 0.9 Mm), with the velocity field overplotted at f = 
243 min and (b): t - 286 min. (c); Fieldlines traced around the center of the EFR at f = 243 min, showing the arcade. Overplotted is 
the full magnetic field vector colored by the component of the field at z = 0.9 Mm. Red is positive and blue is negative. The red 
rectangular inset marks the same region as the one in (b). (d): Fieldlines at f = 286 min, showing the deformation of the PIL. (e): 
Isosurface of > 60 km s ' showing the V-shaped jet at f = 248 min. (f); The emission of several jets at f = 286 min (isosurface 
of Vj > 60 km s“'). 


at the interface. Thus, we find that the deformation of the PIL 
leads to the parallel emission of new jets from the interface. The 
highly intermittent evolution of for t > 280 min in Fig. is 
due to this effect. 

The 3D structure of the emitted jets before and after the de¬ 
formation of the PIL is shown in Fig. I5.IZF- In Fig. |7^ we show 
the 3D structure of the t - 248 min jet, by visualizing the total 


velocity field as an isosurface (red color) with V > 60 km s '. 
At low heights, the reconnection upflow has the shape of a sharp 
wedge with its length oriented vertically, above the reconnec¬ 
tion site. At larger heights, it becomes confined by the envelope 
field and becomes diverted into two lateral flows, following the 
reconnected fieldlines. Thus, the overal configuration of the up¬ 
ward jet adopts a V-like shape. The magnetic fieldlines (white). 
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Fig. 8: Visualization of the heating term (H) at f = 268 min for plasma having temperatures between (a): 0.01 - 1 MK and (b): 
1 - 10 MK. 


which have been traced from within the lateral parts of the jet 
illustrate the envelope field. We also identify the two magnetic 
bipoles and the arcade underneath the reconnection site. At this 
time of the evolution, the downward reconnection jet is much 
weaker and, thus, it appears as a horizontal ribbon-like structure 
just above the arcade. Overplotted is the full magnetic field vec¬ 
tor (horizontal cut) at z = 0.9 Mm. 

In Fig. 1^, we visualize the total velocity field as an isosur¬ 
face with V > 60 km s ' at f = 311 min. The striking difference 
is the appearance of several jets (here, four jets) from the inter¬ 
face. All the upward jets have similar V-like shapes. The parallel 
emission of several hot jets from the interface starts after the de¬ 
formation of the PIL and it is followed by intermittent heating 
of the plasma at the interface. The latter is discussed more in the 
following section. 


3.2.3. Intermittent heating at the solar atmosphere. 

To study how much and where the dense plasma is heated during 
the emission of the jets, we calculate the term; 


H = 



(9) 


We do the calculation for temperatures in the range of 0.01 - 
1 MK and 1-10 MK, as a preliminary attempt to separate the 
cooler (i.e. < 1 MK) from the hotter component of the plasma 
(i.e. > 1 MK). Hereafter, for simplicity, we call this term the 
“heating term”. Figure]^ shows the heating term after the emis¬ 
sion of the first four successive jets (i.e. at f = 268 min), during 
which the temperature enhancement of the associated plasma is 
more pronounced (see Fig.|^). 

In the lower temperature regime (Fig. |^), we find that heat¬ 
ing is more evident in a fan-like arcade, which extends (verti¬ 
cally) from the chromosphere to the low corona. The arcade is 
formed by the reconnected fieldlines, which are released down¬ 
wards, after the onset of the jets at the interface. The fan-like 
shape is formed naturally due to the 3D lateral expansion of the 
interacting magnetic bipoles and it is in agreement with the 3D 
shape of the arcade shown in Fig. |7};. 

In the higher temperature regime (Fig.[^), the heating is not 
smoothly distributed across the arcade, but it is much more in¬ 
termittent. There are four small scale (< 1 Mm) brightenings, 
which all occur within the area between the chromosphere and 
the transition region. These brightenings occur exactly under 
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the reconnection sites of the four bi-directional jets, which were 
emitted during the t - 240 - 260min period. The reconnection, 
which triggers the jets, occurs mainly in the corona. However, 
the downward jet(s) collide with the top of the arcade forming 
a termination shock. There, the plasma is compressed locally at 
high temperatures, causing heating of the dense material. 

A very similar process, regarding the intermittent heating, 
was reported in recent more realistic simulations (radiative trans¬ 
fer, MHD) numerical experiments ( Archontis & Hansteen|2014) l. 
In those simulations, small-scale flaring activity was produced 
at the current sheets between interacting magnetic loops. The 
small flares were triggered by plasmoid-induced reconnection at 
the interface current layers. In the present work, we are witness¬ 
ing an analogous process: patchy reconnection associated with 
the ejections of plasmoids leads to plasma heating, which re¬ 
sembles the onset of nano-micro flares in the model by | Archom] 
tis & Hansteen (2014[). The present model does not include all 


the necessary energy terms (radiative transfer, heat conduction, 
radiation, etc.) for a complete treatment of the thermodynami¬ 
cal behavior of the system. The lack of these terms affects the 
estimated value of temperature (e.g. of the jets) in the present 
experiments. 


3.2.4. Energy transfer and magnetic flux 

To study the transfer of flux into the solar surface and above, we 
first calculate the Poynting flux. We split it into the “emergence 
term”, which is the part of flux representing the direct emergence 
of the magnetic field, and the “shear term”, which represents the 
flux due to the horizontal motions ( |Kusano et al.|200^ |Magara| 
|& Longcope|2003| l: 


1 r 2 2 

^emergence — J (^x 




Fshear — . I i^x^x 4" ByVy^B^dxdy. 


( 10 ) 

( 11 ) 


We calculate these fluxes at the mid-photosphere (zo = 0.8 Mm), 
the base of the transition region (zq = 2 Mm) and the low corona 
(zo = 3.2 Mm). 

In Fig. 1^, we plot the emergence term (solid line) and the 
shearing term (dashed line). At the mid-photosphere (black), we 
find the emergence term to be greater than the shearing term 
between f = 140 - 155 min, as expected during the beginning 
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Fig. 9: (a): Temporal evolution of the Poynting Flux. Solid lines show the emergence term, Femergencej and dashed lines show the 
shearing term, Fshear- (b): The unsigned magnetic flux (c): The axial flux for the bipole’s center in the positive ^-direction 
(bipole AB). (d): Flux difference The dashed vertical line marks the initiation of the undulation of the PIL and the dotted 

vertical lines show the times for the onset of the hight-speed jets. 


of the emergence. The shearing term becomes greater than the 
emergence term, when the inner polarities collide. Owing to this 
collision, plasma is pushed sideways along the PIL and towards 
the T-direction, contributing to the shearing term (e.g. Fig. [^, 
flows along the PIL). Both terms reach their maximum value at 
t =» 200 min and soon thereafter start to decrease. The drop in the 
two terms is due to the fact that the dynamical emergence and ex¬ 
pansion of the field becomes less pronounced at the later stage 
of the evolution of the system. At larger heights (blue lines), we 
find that the shearing term is more dominant due to the expan¬ 
sion of the field. 

To study the vertical flux transfered at the same heights, we 
calculate the unsigned magnetic flux through the Ti^-plane, 

= f IB.Idxdi/. (12) 

At the photosphere, increases monotonically until t - 
255 min (Fig. |^, black line). The drop of the flux thereafter 
is evidence of flux cancellation at the photosphere, even during 
the deformation of the PIL after t - 260 min. Measuring (t>xy at 
larger heights (red and blue lines) shows that less vertical flux is 
transported higher, as expected. 

We also calculate the normalized axial flux through the xz- 
plane at a position yo, from the heights (zo) up to the top of the 
numerical box (Zmax), i-e. 

Bydxdz / I By(t - 0)dxdz. (13) 

zo Jo 

We calculate this axial flux at yo = +6 Mm, which are located 
around the center of the bipoles. We find the same flux at both 


locations, meaning that the same amount of axial flux emerges 
at both bipoles. In Fig.[^, we plot this flux for the bipole center 
in the positive (/-direction. An interesting result is that for each 
bipole, only 12% of the initial axial flux emerges above the pho¬ 
tosphere (black line). From this, only 7.5% gets above the corona 
(blue line). Adding the axial flux of the two bipoles, a sum of 
24% of the initial flux emerges above the photosphere and 15% 
above the corona. Thus, this study shows that non-twisted sub- 
photospheric fields with these characteristics are not capable of 
transfering a considerable amount of flux above the photosphere. 

Owing to the non-twisted magnetic field in our initial condi¬ 
tions, there is only positive By at the beginning of the simulation. 
The emergence and expansion of the field forms the B^ and B^ 
components. In our configuration, negative By component at the 
Tz-midplane could be formed through reconnection. This com¬ 
ponent is part of the arcade’s magnetic field lines. In order to 
estimate how much negative By is created on this plane, we cal¬ 
culate the normalized flux difference: 

^ \ (l^»l “ ^y) ^ 0)dxdz. (14) 

Here AO^^ measures only the flux of negative By through the 
Tz-midplane, from a height zo up to the top of the numerical 
box Zmax, normalized by the initial positive flux before the emer¬ 
gence. So, A<t)^j indicates how much of the initially positive flux 
is transformed into negative flux, through reconnection. This cal¬ 
culation is valid only at the xz-midplane and due to the absence 
of twist. 

Figure]^ shows above the mid-photoshere (black line), flux 
difference increases after t - 170 min. This is due to reconnec¬ 
tion at the interface. This leads to the formation of new fieldlines 
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(a) )v=10, t=3min 
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Fig. 10: First row: z component of the forces along the vertical direction (x — 0, y - 0) at t - 3 min, for (a) d = 10 and (b) d = 5. 
Second row: Temporal evolution of (c) z-component of the magnetic tension and (d) magnetic pressure, measured at the center of 
the flux tubes. We plot the quantities until the flux tubes reach the photosphere. 


that point in the negative [/-direction, joining the two inner polar¬ 
ities, and forming the arcade. We find that increases rapidly 
after the major ejections (dotted vertical lines), peaking after the 
t = 240 - 260 min series of jets. This maximum shows that 2.6% 
of the initial By flux of the system is transformed into negative 
By. In the corona (blue line), we only find 0.3% of the flux. As a 
result, most of the post-reconnection flux at the interface resides 
between the photosphere and the low corona. 

After f = 260 min (vertical dashed line), we find that 
starts to decrease drastically. This is when the deformation of 
the PIL starts to occur (Fig. |^, ^). At this moment, the inter¬ 
face cannot be studied by using the 2D vertical .rz-midplane, and 
A<[);cj does not represent an accurate measurement of the nor¬ 
malized flux difference. This happens because, as we described 
previously, the field lines of the arcade follow the deformation 
of the PIL (Fig. |7}i). Thus, some of the field lines which were 
previously perpendicular to the interface, they are now pointing 
towards the x-direction due to the undulation (e.g. Fig.jTJ), lines 
around x - +0.3 Mm would tend to be parallel to the +-axis). 
As a result, their By (B^) component decreases (increases), and 
consequently, AO;^^ decreases. Also, the field lines that do not 
rotate, tend to move away from the .rz-plane due to the flows in 
the [/-direction (e.g. around the EFR center, PIL position is lo¬ 
cated + = 0, [/ 1.3 Mm). This further decreases By on the 

TZ-midplane, and contributes to the decrease in AO^^j.. 


4. Case of /I = 5 and comparison with /I = 10 

In the following, we describe the emergence of the non-twisted 
flux tube with T = 5, and make a comparison to the T = 10 case. 
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4. 1. Emergence from the solar interior to the photosphere 

In the following, we firstly show the differences between the two 
A cases when both flux tubes are located in the solar interior. In 
our experiments we apply a density deficit perturbation, which is 
practically controlled by A (Eq.0. In order to study how differ¬ 
ent A affects the initial emergence from the solar interior to the 
photosphere, we examine the cross section of the flux tubes at 
the vertical Tz-midplane (i.e. y = 0). There, both flux tubes are 
equally buoyant (same Ap in Eq. since the pressure balance 
terms are the same for both cases). 

At f = 0 min, the two flux tubes have equal forces applying 
to them (at y - 0) along the vertical z-axis. In Eig. [T0^,b we 
plot for /I = 10 and 4 = 5 respectively the z-component of the 
forces along the z-axis at f = 3 min. We find that: (i) the gas pres¬ 
sure force minus the gravitational force (green line) and (ii) the 
magnetic pressure force (orange line) are approximately equal 
in both cases. On the other hand, the magnetic tension force (red 
line) shows a distinctive difference from the beginning of the 
emergence. Eor 4 = 5, the tension is higher because the apex 
of the buoyant part of the tube is more convex. Thus, the lower 
tension case experiences higher total force (black), making the 
flux tube to reach the photosphere faster than the one with higher 
tension and lower total force. Indeed, by comparing the apex po¬ 
sitions of the flux tubes, we find that the 4=10 flux tube reaches 
the photosphere at f = 26 min while the T = 5 at f = 34 min. 
Eig.p^ shows the temporal evolution of the z-component of the 
magnetic tension at the center of the flux tubes, until the flux 
tubes reach the photosphere. We find that the tension is up to 
3 times higher for smaller A (dashed line) throughout the initial 
emergence. As a result, the high tension of smaller A continu¬ 
ously decelerates the emergence of the flux tube. 
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Fig. 11: Forces and flux tube cross sections of d = 10 and d = 5, before they reach the photosphere (at t -26 min and at f = 34 min 
respectively). (a,d): Pmag of the flux tubes at the y - 0 cross section. The solid white lines (horizontal and vertical) pass through 
the flux tubes centers. Black P^ag contours have the same values in both panels. Vertical white dashed lines mark the position of 
the flux tube’s outer part. (b,e): z-component of forces along the vertical direction through the flux tube’s center (direction marked 
with the vertical solid white line in (a,d)). Green is gas pressure force minus gravity force, orange is magnetic pressure force, red 
is magnetic tension force and black the total force. (c,f); The x-component of forces along the horizontal direction through the flux 
tube’s center (direction marked with the horizontal solid white line in (a,d)). 


This faster emergence of the higher A flux tube, causes it 
to expand more in the vertical direction. Therefore, it leads to a 
faster drop of its gas and magnetic pressure in comparison to that 
of d = 5. In Fig. 101, we plot the magnetic pressure P^ at the 
center of the flux tubes. We plot this quantity until the moment 
each flux tube reaches the photosphere. Indeed we And that for 
d = 5 (dashed line), the flux tube has higher magnetic pressure 
at its center than for d = 10 (solid line). 

Figures [TT^ and Fig. [TT] l show the magnetic pressure at the 
xz-midplane (for d = 10 and d = 5, respectively) when the flux 
tubes reach the photosphere (at different times). The solid white 
lines indicate the vertical and horizontal directions through the 
flux rope’s center. As we showed before, the faster vertical emer¬ 
gence leads to lower internal pressure for the higher d flux tube. 
This is apparent again by examining the z and x components of 
the forces along the vertical and horizontal direction through the 
center of the flux tube (Fig. [TTJ),c for d = 10 and Fig. [TT^,f for 


d = 5). 

Along the vertical direction, the comparison of the forces 
(Fig.[n]r, e) shows again that the magnetic tension is higher for 
d = 5. We find that both the gradient of gas (green line) and 
magnetic pressure (orange line) are lower for d = 10, indicating 
the faster drop of pressure due to the faster emergence. 

The comparison of the horizontal expansion of the flux tubes 
(dashed lines in Fig. [U ,d) indicates that d = 10 expands more. 
Indeed, the x-component of the total force in the horizontal di¬ 
rection (black line Fig. [U ,f) shows that total force is higher 
for d = 10. The forces that mainly act on the flux tube in the 


horizontal direction are the outward magnetic pressure, trying to 
expand the flux tube, and the inward gas pressure. In our initial 
flux tube, we have no twist, so initially there is no tension in the 
transverse direction. Still some inward tension is developed (red 
lines) due to the expansion of the flux tube and possibly due to 
internal reconnection of the held lines during this initial emer¬ 
gence. 

An important result arises from the interior of the flux tubes 
during the initial emergence. Notice the contours in Fig. [TT}i,d. 
It is apparent that the d = 5 case emerges coherently. On the 
other hand, the d = 10 case develops two minor lobes in its in¬ 
terior. The lobes are formed as the center of the emerging flux 
rope starts to split in two parts. This also affects the plasma 
density in these regions. Density becomes lower in the vicin¬ 
ity of the lobes (due to higher magnetic pressure) and higher in 
the center of the flux rope (where the magnetic pressure drops). 
As a result we And outward directed gas pressure force (i.e. to¬ 
wards the lobes) in the interior of the d = 10 flux tube (around 
-0.4 Mm<x<0.4 Mm, Fig. my 

A similar behavior has been reported in previous studies of 
emerging flux tubes (e.g.|Schuessler|1979[ Longcope et al. 1996[ 


[Moreno-Insertis & Emonet|1996| l. Those studies showed that the 

cross section of rising flux tube was distorted through the in¬ 
teraction with the background non-magnetized medium. More 
precisely, they found that the emerging flux tube displaced the 
plasma above it and forced it to flow from the apex to the flanks 
of the tube. This created a drag that deformed the cross section 
of the flux tube and eventually created two counter-rotating vor- 
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Fig. 12: Cross sections of P^ag at the xz-midplane for (a): d = 5, (b): A = 10, (c): A - 15, (d): A - 20. Purple is zero and dark red 
is 3.6 X 10"^ erg cm“^. Black lines are P„,ag contours (same contour values for all images). Gray contours are the [/-component of 
vorticity (same contour values for all images). Solid lines are positive values and dashed lines are negative values. 


tices. In our numerical simulations (e.g. d = 10 case) the ap¬ 
pearance of the two minor lobes just before the flux rope reaches 
the photosphere (Fig. [TT^). On the other hand, it seems that in 
the d = 5 case this lobe formation is inhibited and the flux rope 
emerges as a whole, without distortion (Fig.[TT]l). 

To further investigate the possible role of d in the coherence 
of the flux tube during the emergence in the solar interior, we 
also performed simulations for d = 15 and d = 20. In Fig. 12 we 


plot the magnetic pressure of the cross section of the flux tubes 
on the xz-midplane before the flux tubes reach the photosphere, 
for each d. The black contours show P,„acj- To do a closer com¬ 
parison between the different d cases, we show the same value of 
Pmag- The coutours of Pmag indicate that lower d results in more 
coherent flux tubes that carry higher magnetic pressure flux at 
the photosphere (see dark red color inside the flux tubes). Also, 
we And that flux tubes with higher d are more susceptible to split¬ 
ting into two lobes. Indeed, the d = 20 case (Fig. [T^) develops 
two very distinctive lobes, while in the d = 10 (Fig. [T^) those 
lobes are marginally separated (the dark red region is more elon¬ 
gated rather than separated). This splitting is caused (in a similar 
manner to the afore-mentioned studies) by plasma flow from the 
apex to the flanks. This flow creates an aerodynamic gas pres¬ 
sure difference below the regions where we And the maximum 
horizontal expansion (around z = -1 Mm). Then, the pressure 
difference drives plasma motions. To illustrate this, we overplot 
the contours of the [/-component of vorticity (gray contours in 
Fig. 12 1 . We indeed And two counter rotating plasma flows along 
the flanks of the flux tubes. We note that the same vorticity val¬ 
ues (for instance the second contour from inside out) are located 
deeper into the flux rope for larger d. This signifies the presence 
of stronger gas-pressure driven plasma flows, which eventually 
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deform the flux tubes. We And that lower tension along the flux 
tube axis (higher d) results in more pronounced splitting of the 
flux tube. We should also mention (not shown here) that the flux 
tubes with d = 15 and d = 20 failed to emerge above the photo¬ 
sphere during the running time of our simulations. 

We now focus on the d = 5,10 cases, which eventually 
emerge above the photosphere. We have found that the higher 
tension flux tube reaches the photosphere later and without lower 
expansion. Consequently, the flux tube with higher tension reach 
the photosphere with higher magnetic pressure (i.e. higher mag¬ 
nitude of magnetic held) and lower plasma This has very im¬ 
portant consequences on flux emergence at and above the solar 
surface. We And that the d = 5 flux tube satisfies the instability 
criterion (Eq. earlier due to its higher magnitude of magnetic 
held and the lower fi. Since the emerging held is non-twisted, 
it spreads horizontally and it becomes buoyantly unstable at two 
regions (similar to the d = 10 case), forming two bipoles. A main 
difference compared to the d = 10 case is that these bipoles ap¬ 
pear at the photosphere earlier, (at t - 66 min, f = 154 min for 
d = 10) and at a shorter distance (at y - +1.8 Mm, y - 4.9 Mm 
for d = 10). 

In the following sections we are going to highlight how the 
higher tension tube eventually goes through a more dynamic 
emergence and transfers more Poynting and axial flux above the 
photosphere than the lower tension tube. 


4.2. Jet and associated plasma dynamics 

The dynamics of the jets depends on the magnetized plasma 
properties at the interface. We And that for d = 5, the gradi- 
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Fig. 13; (a): Temporal evolution of maximum (black) and minimum (red) I 4 above the photosphere. The vertical dashed lines 
correspond to the plasma ejections at f = 114, 120, 132 min. (b): Maximum temperature and density around the locations of 


maximum V^. 



Fig. 14: (a): Fieldline topology at the vertical j/z-midplane, showing the plasmoid(s) at the interface (at t - 108 min). The horizontal 
cut shows the component of the magnetic field at the photosphere, and the purple-green colorscale shows density in the yz- 
midplane. The fieldlines are colored according to the values of By. (b): Isosurface of V > 60 km s ' showing the f = 114 min jet. 
White fieldlines show the stretch of the envelope field by the ejection. Overplotted is the full magnetic vector colored by the B^ 
component of the field at z = 0.9 Mm. 


ent of the magnetic field across the interface is much higher. The 
resulting curi'ent sheet is 3-5 times stronger. This leads to dif¬ 
ferent characteristics (e.g. higher plasma temperatures) for the 
produced plasma jets. 

Fig. [T3^ shows the emission of three high-speed (200 - 
600 km s^) bi-directional jets from the interface, followed by 
lower speed jets (< 100 km s“' after t « 135 min). The absolute 
values of the speeds of these jets are, on average, higher than 
in the T = 10 case. This is due to the higher field strength (and 
lower plasma yS) of the magnetized material within and around 
the interface. However, it is worthwhile mentioning that in both 
cases the speeds of the reconnection jets are comparable to the 
local Alfven speed, as expected from the theory of reconnection. 
In this manner, we find that the speeds of the jets in the two cases 
are similar (determined by the Alfven speed). 

Regarding the heating due to reconnection-driven jets, the 
plasma of the upward jets has a higher temperature in the 4 = 5 
case. This is shown in Fig. [T^, which illustrates the temporal 
evolution of the maximum temperature at the sites of maximum 
Vj. The temperature can reach values up to 15 MK (first jet). The 
following jets have progressively lower temperatures. In general. 


the temperature increases by a factor of 1 ![}. As we discussed 
earlier, plasma-/? around the reconnection site is higher when A 
is smaller. In the time period after the onset of the jets (i.e. f > 
240 min for 4 = 10, f > 105 min for 4 = 5), we calculated that 
(on average) /? a; 0.015 for 4 = 5 and =» 0.17 for 4 = 5. This 
explains the differences in the temperature of the jets when 4 is 
varied. 

The corresponding density evolution shows a very good cor¬ 
relation with the temperature variations (maxima of temperature 
and minima of density occur at the same time). In comparison 
to the 4=10 case, the density variations lie (on average) in the 
same regime 10“*® - 10“'"* g cm“^. Thus, in both cases, the jets 
bring dense plasma into the solar atmosphere. 

4.3. Fieldline topology and 3D structure of the jets. 

In both cases, the process that leads to the onset of the jets is sim¬ 
ilar. Namely, the interaction of the two bipoles leads to the for¬ 
mation of a strong current layer at the interface. The current layer 
becomes subject to the tearing-mode instability and plasmoid- 
induced reconnection leads to the emission of hot and high-speed 
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jets. Hereafter, to avoid repetition, we describe briefly only the 
topology during the first jet. 

Figure 


14 


shows the fieldline topology at f = 108 min, a 
few minutes before the emission of the first jet. The vertical 
midplane shows plasma density (blue-red colorscale), and the 
horizontal slice is at the photosphere (black-white). We note 
that now the fieldlines are colored by the By component of the 
magnetic field. In our non-twisted configuration, the succession 
of positive and negative By highlights the locations of the plas- 
moid(s) and the arcade. At this time, one plasmoid has already 
been ejected upwards and another one has been formed at lower 
heights within the interface. The ejected plasmoid has collided 
with the envelope field and, consequently, it undergoes severe 
deformation. The ejection of the plasmoid is confined by the 
downward tension of the envelope fieldlines. The other plasmoid 
is ejected soon thereafter and it is also trapped by the closed am¬ 
bient fieldlines. 

After the ejection of the plasmoids, several jets are emitted 
from the interface. Fig. [T4] r shows the 3D structure of the first 
jet at f = 114 min. The upward reconnection jet adopts a V- 
like shape (similar to the 4=10 case). During its emission, the 
fast jet stretches the ambient field vertically but it does not break 
out. The jet is diverted into two lateral flows, which are running 
along the reconnected fieldlines of the envelope field. This leads 
to an apparent heating of the envelope fieldlines that join the 
outer polarities of the system. 

Just underneath the upward V-like jet, there is another fast 
jet, which is less extensive and it adopts an inverse V-like config¬ 
uration. It extends downwards until the top of the arcade, which 
is formed by the O-like fieldlines that are rooted to the minor po¬ 
larities of the system. The jet compresses the plasma at the top of 
the arcade and heats it to « 15 MK. The further evolution of the 
system is qualitatively similar to the 4=10 case. For instance, 
the emission of several plasmoids from the interface leads to the 
onset of a series of jets, which can exist in parallel (e.g. like in 
Fig-IZf)- Also, we find that the PIL undergoes deformation, de¬ 
veloping horizontal undulations, around the center of the EFR, 
after f = 119 min. 


also find evidence of flux cancellation at the photosphere (after 
t — 120 min). We find a less profound decrease of flux above 
the transition region, while it is not clearly apparent above the 
corona. In Fig. [TS] ; we plot the normalized axial flux of one 
bipole (the bipole that is located in the positive j/-direction) 
above different atmospheric heights (solid lines). For compar¬ 
ison, we plot with dashed lines the 4 = 10 values above the 
same heights. Above the mid-photosphere, we find that 26% of 
the initial axial flux emerges and from it 16% is located above 
the corona. Adding the axial flux of the two bipoles, we find 
that 52% of the initial axial flux emerges above the photosphere, 
which is 1.5 times the 4 = 10 axial flux. Above the corona, we 
find 32% of the axial flux to emerge, which is 2.4 times larger 
than in the 4=10 case. 

These differences are caused by the emergence of stronger 
magnetic field in the lower 4 (higher tension) case. This results 
in transfer of more flux and energy into the solar atmosphere. 
However, we note that in both cases 2/3 of the initial axial flux 
remain below the corona and only half of the axial flux emerges 
above the photosphere. Therefore, we find that in the cases under 
study a large amount of axial flux either stays at the photosphere 
or even submerges. 

We also calculate the flux difference AO^^j.. Above the pho¬ 
tosphere (Fig. 151, black line), AO^^j. increases rapidly after the 


first ejection (dotted vertical line). Its maximum, at f = 120 min 
stands for 5% of the initial flux being transformed into negative 
By. This flux difference contains both the arcade and the plas- 
moid’s By. During the upward motion of the plasmoid(s), some 
of the plasmoid’s lines reconnect with envelope fieldlines. Part 
of this flux is added to the arcade, contributing positive to AO^.^, 
and part of it enhances the envelope field, contributing negative 
to We find that overall, decreases as the plasmoid(s) 
moves upwards. The vertical dashed line indicates the start of 
the deformation of the PIL. Again, as in 4 = 10, AO;^^ decreases 
as the arcade follows the PIL undulation. As we mentioned in 


3.2.4 Afi) indicates the transformation of flux through re¬ 


connection. Thus, we find increased reconnection along the lobe 
interface in the smaller 4 case. 


4.4. Energy Transfer and Magnetic Flux 

In Fig.ffS^, we plot the shearing (dashed) and emergence (solid) 
term of Poynting flux, above the mid-photosphere (black), tran¬ 
sition region (red) and the corona (blue). Energy is injected 
above the photosphere after f = 60 min, much sooner than the 
t = 140 min of 4 = 10, due to the faster emergence of the flux 
tube above the photosphere. One difference between the two case 
is that more energy is injected into the atmosphere in the smaller 
4 case. Indeed, the emergence term in 4 = 5 is 4.2 times higher 
that the one in 4 = 10, while the shearing term is 2.4 times 
higher. This is because the emerging field with larger tension 
(i.e. 4=5 case) brings stronger magnetic field at the photosphere 
and evolves/expands more dynamically (e.g. faster) in the so¬ 
lar atmosphere, as we showed in the previous sections. As a re¬ 
sult, more flux/energy is transfered into the outer atmospheric 
layer, when the subphotospheric emerging field has larger ten¬ 
sion. Therefore, our experiments reveal that this specific param¬ 
eter of the subphotospheric emerging field is important towards 
understanding the dynamics of the emergence at and above the 
solar surface. 


The unsigned flux, (^xy (Fig- 15 t), indicates that more flux is 
injected into when 4 is small. Ihe maximum of the flux above 
the mid-photosphere (black), transition region (red) and corona 
(red) is 1.8 times larger the corresponding flux of 4 = 10 . We 
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5. Summary and Discussion 

In the present work we studied the emergence of a non-twisted 
flux tube in a highly stratified atmosphere. The parameter we 
varied in our experiments was the buoyant part of the flux tube 
(4). 

We examined the initial emergence stage (below the pho¬ 
tosphere). We found that the magnetic tension of the emerging 
part of the rising flux tube (smaller tension for higher 4) has a 
significant effect on the emergence process. More precisely, the 
4=10 tube is distorted into two smaller lobes. Eor 4 = 5, the 
tube emerges in a more coherent manner. Moreover, smaller ten¬ 
sion leads to higher upwards total force. As a result, in the 4=10 
case, the tube reaches the photosphere earlier and expands more 
horizontally, while it is still within the subphotospheric region. 
This faster emergence/expansion drops the flux tube’s internal 
magnetic pressure faster 4 = 10. Subsequently, the higher 4 
(smaller tension) case reaches the photosphere having smaller 
magnitude of the magnetic field. Eventually, this difference in 
the magnitude of the magnetic field, makes the smaller 4 (higher 
|B|) to emerge above the photosphere faster, although it reaches 
the photosphere later than the higher 4 case. In both cases, the 
initial flux tube becomes buoyant in two positions along the 
tube’s length, forming two aligned bipoles at the photosphere. 
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Fig. 15: (a): Temporal evolution of the Poynting Flux. Solid lines show the emergence term, Fg^„gc„ce, and dashed lines show the 
shearing term, Fshear- (b): The unsigned magnetic flux (c): The axial flux for the bipole’s center in the positive i/-direction 
(bipole AB) (solid). Dashed lines stand for the axial flux of T = 10. (d): Flux difference The dashed vertical line indicates 

then beginning of the undulation of the PIL and the dotted vertical lines indicate the times of the onset of high-speed jets from the 
interface. 


Since the d = 5 emerges first, the horizontal distance between 
the two bipoles at the photosphere is smaller. 

Simulations with A - 15,20 indicated that distortion be¬ 
comes more intense for higher A. The origin of this distortion 


seems to be similar to the ones described by e.g. Schuessler 
(1979i; [Longcope et al. ( |1996| l; [Moreno-Insertis & Emonei 


(1996i. Flows at the flanks of the emerging flux tube increase 
the plasma vorticity locally. It is likely that the drag forces at the 
flanks of the flux tube cause the distortion. Although the results 
presented here are preliminary, we do find indications that the 
convex apex of the emerging flux tube helps the magnetic field 
to rise coherently. A full analysis with different magnetic field 
strengths, initial depths of the flux tubes and various twists is 
necessary for further conclusions. 

The magnetic bipoles interact and reconnect at the interface 
between them, forming an envelope field and an arcade. We 
find that the interface current sheet is « 5 times higher when 
d = 5 which leads to more effective (faster) reconnection. From 
the interface current sheet, cool plasmoids are initially ejected, 
followed by several reconnection jets which are running with 
speeds comparable to the local Alfven speed. All upward jets 
are confined by the envelope field and adopt a V-like shape. The 
downward jets compress the plasma at the top of the arcade and 
heat it to a few MK. The heating occurs in an intermittent manner 
due to the parallel emission of the jets at the interface. A simi¬ 
lar process has been reported in more realistic experiments re¬ 


garding the onset of nano-micro flares by Archontis & Hansteen 

( |MT4t . 

We have also found that the value of d has implications on 
the energy and flux input into the atmosphere. Indeed, the emer¬ 


gence and shearing term of the Poynting flux are 4.2 and 2.4 
times higher in d = 5. The former is due to the faster emergence 
and the latter is due to more pronounced horizontal flows (along 
the PIL) at the interface. The vertical flux transfer is 1.8 times 
larger in d = 5. In d = 5 (d = 10), 52% (30%) of the initial 
axial flux is transfered above the photosphere, and 24% (15%) 
above the low corona. The increased energy and flux transfer for 
smaller d occurs as in this case stronger magnetic field is brought 
at the photosphere. This causes an overall faster and more dy¬ 
namical evolution in the solar atmosphere. 

In our aligned bipole configuration, we noticed that after the 
plasmoid ejections the PIL starts to gradually deform. The erup¬ 
tions induce a total pressure deficit (at various locations) along 
the PIL. The flows towards the low-pressure regimes deform 
the PIL, forming finger-like features that are mainly oriented 
along the (/-axis (transverse to the PIL). This makes the post¬ 
reconnection arcade to become a highly curved structure con¬ 
sisting of twisted fieldlines. 

Previous studies that compared cases with different d (e.g. 
[Manchester et al.||2004[ [Archontis & Hood||2010] l showed that 
the buoyant part of a twisted flux tube affects the eruption of 
newly formed flux ropes, which are triggered by shearing and re¬ 
connection within the emerging flux region. In the present study, 
we investigated the effect of d on the emergence of a non-twisted 
field. Although the initial magnetic field configuration and initial 
conditions in our experiments are idealized, the configuration of 
the emerging field at the photosphere and the subsequent dynam¬ 
ical evolution become very complex. The emission of several 
jets from the same interface and the intermittent heating of the 
plasma at different atmospheric heights is just one example of 
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the arisen complexity. Therefore, this parametric study showed 
that the variation of only one parameter of the system can affect 
the dynamics of the system considerably. More work is needed 
to investigate the dynamics of the jets and the atmospheric re¬ 
sponse to the emergence of magnetic fields from the interior of 
the Sun. 
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